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The nuclear shape correction to the g factor of a bound electron in IS-state is calculated for a 
number of nuclei in the range of charge numbers from Z — 6 up to Z — 92. The leading relativistic 
deformation correction has been derived analytically and also its influence on one-loop quantum 
electrodynamic terms has been evaluated. We show the leading corrections to become significant 
for mid-Z ions and for very heavy elements to even reach the 10~® level. 

PACS numbers: 



The ever-increasing precision of measurements and 
theory of the g factor of a bound electron has recently de- 
livered a new value for the electron mass [TJ [5] , and keeps 
providing stringent tests for quantum electrodynamics 
(QED) in strong fields [THS]- It also allows to access 
electromagnetic properties of nuclei such as charge radii, 
as demonstrated in a very recent proof-of-the-principle 
study with a Si^^"*" ion [3], or, as suggested theoretically, 
magnetic moments [1]. Also, it is anticipated that g fac- 
tor studies will yield a value for the fine-structure con- 
stant a that is more accurate than the presently estab- 
lished one when extending the experiments to elements 
with a high charge number Z [5]. 

In a few years, measurements with the heaviest ele- 
ments will be possible [5j. As higher-order nuclear and 
QED contributions to the theoretical value of the g fac- 
tor are strongly boosted with increasing Z, at the present 
10~^° level of relative experimental accuracy [3j[7] or even 
below, our present understanding of atomic structure will 
not be satisfactory. In such strong Coulomb fields, nu- 
clear effects beyond a simple spherical model arise. Fur- 
thermore, QED and nuclear structural contributions are 
intertwined. 

In this Letter we consider the nuclear shape effect, and 
find that while it can be safely neglected in predictions 
for \ow-Z systems, it greatly influences the g factor value 
already for mid-Z elements. At high nuclear charges, its 
inclusion in the theoretical description is mandatory; for 
example, for U^^"*", its relative contribution to the total 
g factor reaches the 10~^ level. We furthermore evaluate 
mixed nuclear-QED terms, i.e. the nuclear shape effect 
on the one-loop QED terms of self-energy (SE) and vac- 
uum polarization (VP). Even these contributions will be 
highly relevant for the interpretation of experimental val- 
ues to be obtained within a few years [5]. Furthermore, 
a comparison of theory and experiment may even yield 
more accurate values for nuclear shape parameters, rel- 
evant in explaining shape phase transitions in nuclear 
structure theory |8 . 

We account for the nuclear quadrupole and hexade- 
capole deformation and derive a formula describing the 
nuclear shape correction to the g factor of a bound elec- 



tron in hydrogen-like ions in the IS'-state. Then we focus 
on systems where the nucleus is spinless and in its ground 
state. Let us start with the definition of the electron g 
factor (we use units with c ^ 1, h = 1, a — e^/(47r), and 
with e = — I e| unless otherwise stated) 
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where 6E stands for the energy correction due to the 
coupling of the electron spin operator s to the external 
magnetic field B, m is the mass of the electron, and (• ■ •) 
stands for the expectation value. The nuclear deforma- 
tion (ND) correction to the g factor, which we present, is 
closely related to the finite-size (FS) effect, which is well 
known in the context of atomic levels. The correction 
emerges when one accounts for a deviation of a nucleus 
from the spherical shape. 

Let us consider now a relativistic Hamiltonian of the 
bound system of an electron and a nucleus in the presence 
of an external, constant and homogeneous magnetic field. 
The interaction of the nucleus with the external magnetic 
field is negligible so that we can write 



H = Ha + Hn + 5Vc - eaA , 



(2) 



where A — ^{B x f), Ha is the atomic Dirac-Coulomb 
Hamiltonian Ha — dp + /3m — r is the vector describ- 
ing the electron's position with respect to the mass-center 
of the system, H^ the (unknown) nuclear Hamiltonian 
and SVc the deviation from the monopole term expansion 
of the exact Coulomb potential of our system, i.e. 
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Now we formulate the framework for the perturbative 
calculation of the energy correction due to the presence 
of the magnetic field and the deviation from the Coulomb 
potential. We start with the equation 

{Ha+HN+5Vc-edA) \^>) = {E^^^+E^°^+5E) j*) , (4) 



where | Vf) denotes the perturbed state of the system, Ea 
the unperturbed ground-state energy of the atom, E 
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the unperturbed ground-state energy of the nucleus and 
5E denotes the energy correction arising due to 5Vc and 
the interaction with the external magnetic field. Next, 
we perform the following approximation: we assume that 
the potential 5Vc does not induce any nuclear transition, 
which means that we neglect the nuclear polarizability 
(NP) effects. However, the derivation of the NP cor- 
rection as well as its numerical values for a number of 
elements can be found in Ref. [S]. In principle there are 
also mixed terms containing both: the NP and nuclear 
deformation effects, but since these effects are both al- 
ready small we neglect the mixed terms. Then we can 
write l^*) = 10') (g) 1-0), where \i>) stands for the ground 
state of the nucleus, denotes the perturbed atomic 
state 10') = 10) + 1 50) and |i50) stands for the correction 
to the unperturbed atomic state |0). We substitute the 
state 1 4") into Eq. ^ and act on the resulting equation 
with (i/jj. Following an obvious reduction we obtain 



{Ha + 5V{r) - eaA) |0') = + 6E) |0') 



(5) 



where 5V(r) = (-01 <5Vc IV')- 

Now we will express SV{r) in terms of a nuclear charge 
distribution p{r). We do not focus on details here since 
the detailed derivation of the formula ^ can be found 
in Ref. [10]. If the nucleus is in the state of vanishing 
angular momentum / = 0, then the potential 6V{'r) can 
be written in terms of the simple formula 
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where r> = max(r, r'). In Eq. (|6| we have introduced the 
radial charge density, i.e. p{r) = ^ dnpir). The terms 
on the right-hand side (RHS) of Eq. ^ corresponds to 
the respective terms on the RHS of Eq. (|3| . Let us note 
that the resulting potential 5V{r) depends solely on the 
radial variable r. 

In our calculation we consider axially symmetric nuclei 
and employ the two-parameter Fermi charge distribution 
with quadrupole and hexadecapole deformation, i.e. 



p{r) 



N 



1 -I- exp [{r — c)/a] 
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where A'^ is a normalization constant, the half-density 
radius is given in terms of spherical harmonics as c = 
Co(l -t- /32^20 + PiYio) with j32 and being the octupole 
and hexadecapole deformation parameters, respectively. 

The g factor corrections of interest to us for the systems 
described by Eq. ([simight be calculated with the help of 
the results of Ref. [11]. We will use Eq. (27) of Ref. [H], 
which for a IS'-state reads 



4(27 + 1) Evs 
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This formula describes a correction to the g factor due 
to the FS effect for arbitrary radial distributions of the 



nuclear charge. E-ps denotes the energy shift due to the 
FS effect described by the perturbing scalar potential, 
and 7 = \/l — (Za)^. Eq. (|8|) takes into account the in- 
teraction with the external magnetic field in the first or- 
der of perturbation theory, whereas the interaction with 
SV{r) is treated to all orders. In practice, the accuracy 
of Eq. ^ is hmited by the accuracy of Eps ■ In order to 
obtain Eps we use Eq. (17) of Ref. [12]. Namely, for the 
IS'-state 

Efs = [l + {aZffiaZ)] (2aZRmj m, (9) 

where /(aZ) = 1.380-0.162(aZ) + 1.162(aZ)2, and i? is 
the effective radius of the homogeneously charged sphere 
that gives the same energy correction as the original nu- 
clear shape. For states with j = 1/2 
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, (10) 



where (r") = J drr'^ r"p{r)/{Z\e\). In our calculation 
<^r^) is computed by numerical integration, whereas <^r^)) 
is found in the literature [T5] . 

It might be instructive to obtain the approximate for- 
mula expressing R in terms of the nuclear parameters. 
For this purpose, we substitute relations (24)-(26) from 
Ref. [To] into Eq. ( 10 ) and expand the result. In this way 
we obtain 
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The values of R obtained with the help of numerical inte- 
gration are in good agreement with those obtained with 



formula ( 11 ). Finally, the nuclear deformation correction 



to the g factor is defined as 



ND 



<55FS {Z, a, /32, Pa) Sg^s {Z, ^fir^, a, 0, O) , 

(12) 

where Sgps depends on the nuclear parameters via R. 

We also take into account nuclear deformation correc- 
tion to the QED corrections. The QED corrections to 
the g factor of first order in a consist of self-energy (SE) 
and vacuum-polarization (VP) corrections. We mostly 
use the calculation scheme as in [TH [TF] . The SE term 
represents the sum of irreducible, reducible and vertex 
parts: SgsE = Sg-nr + Sg^cd + Sgvcr- The irreducible part 
is given by [IB] 



1 m^{E^) - l°Sm)\n){n\[a x r] 



Ech — En 



(13) 



3 



Here, is an angular momentum of the state |(/)) and En 
the energy of state |n), [a x r]^ describes the interaction 
with the external magnetic field, and Sm is a mass coun- 
terterm. S(-E) denotes the unrenormalizcd self-energy 
operator defined as 

i f (an\I{Lu)\nb) 
2^ J E-oj-En{l-iO)' 
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where /(cj, xi, X2) = e^a^a'^D|^^^{uJ,Xl,X2) with the 
Dirac matrices = and the photon propaga- 

tor D^^ijjj, xi,X2)- The expressions for the reducible and 
vertex parts read |16| 

Sg.cd - —(01 [a X r], 10) (01 -4 S(i?) I 
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(16) 
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(i;^ - w - i;„, (1 - iO))iE^ -oj- En,{i - io)) ■ 

Both the reducible and vertex parts are ultraviolet- 
divergent, whereas the sum Sg^r — Sg^cd + Sgvcr is finite. 
To separate the ultraviolet divergencies, the expression 
(13) is decomposed into zero-, one-, and many-potential 



terms, and the expression (15 1 is decomposed into zero- 
and many-potential terms. All zero- and one-potential 
terms are evaluated in momentum space. The remaining 
many-potential parts Sg^'^j^'^ and 5gil^^ are calculated in 
coordinate space [T7]. The angular integration and the 
summation over intermediate angular projections in the 
many-potential terms Sg^^j^'' and Sgil^'^ are carried out in 
a standard algebraic manner. The many-potential terms 
involve infinite summation over the quantum numbers 
K = ±(j -f 1/2). This summation is terminated at a 
maximum value \k\ = 15 — 20, and the residual part of 
the sum is evaluated by a least-square inverse-polynomial 
fitting. The summation over the Dirac spectrum at fixed 
intermediate k was carried out with the dual-kinetic-basis 
(DKB) method [TS]. 

The VP correction consists of two parts, originating 
from the vacuum-polarization diagram with the magnetic 
interaction inserted into the external electron line (the 
electric-loop contribution Sg^p) and into the vacuum- 
polarization loop (the magnetic- loop contribution 6 gyp). 
We took into account only the leading term since it guar- 
antees the accuracy needed in the present work. The 
leading term originates from the electric loop and can be 
described by the Uehling potential. The VP correction 
is calculated using the DKB method as well. 

To calculate the nuclear deformation correction to 
QED corrections (QED-ND), we used an effective ra- 
dius of the nucleus. We employed the shell, homoge- 
neously charged sphere and Fermi nuclear distributions. 
Although the QED-ND-calculation results are numeri- 
cally stable and do not depend of the nuclear model, they 
should be considered as an estimate. 



TABLE I: The nuclear deformation correction to the g fac- 
tor (5(;nd). The values of {r^Y''^ and their uncertainties 
originate from Ref. [13]. Unless otherwise stated the value 
of a was obtained by equation a = f/(41og3) (Ref. |19]). 
where t — 2.30 fm and rather conservative error bars were 
assumed. T he parameters j32 were estimated with the help 
of formula (17 1 unless a reference is given. For ^^''U and 
^^*U we employed, in accordance with |20| . /34 = 0.08(5) and 
Pi — 0.07(5), respectively; for Si, in accordance with [21j . 
/34 — 0.08(5); for other elements we assumed j34 — 0.0(1). 
The values of cq are set by the requirement that the numer- 
ically computed {r^)^^'^ should be equal to the established 
ones [13]. See also the comments in the text on the uncer- 
tainty of Sqnd ■ 
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-1.28(28) 


10" 


-6 



"Ref. [22]. 
''Ref. Ip. 

"The sign originates from Ref. |21| . 



Table (|T| presents numerical values of i5(7nd obtained 
with formula ( 12 1 for various hydrogenlike ions. De- 
pending on the availability of experimental data, the 
quadrupole deformation parameter /32 was either taken 
from literature or estimated with the help of the for- 
mula [2S1 
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where: _B(i?2;0+ 2^) stands for the reduced prob- 
ability of E2 transition from the ground state 0+ to a 
state 2^, Rs is the effective radius evaluated with the 

help of the simplified formula Rs — \J | {r'^), and the 

values of B{E2; 0+ ^ 2+) were evaluated with the help 
of the data from Ref. [22] . The uncertainty of Sg^-o is es- 
timated as the quadratic sum of uncertainties related to 
the nuclear parameters. We assume rather conservative 
error bars for the parameters a, ^2 and to account for 
an uncertainty caused by the fact that for some nuclei 
the parameters were compiled from various experimental 
results. 

As for the QED-ND effect, it is significantly smaller 
than the leading ND correction. For example, we ob- 
tained the following values: 5(3) -10-9, 4(2)- 10-^ 1.2(4)- 
10"^°, 1.3(1. 1)-10-", 1.5(3)-10-", and 1.6(6)-10"i2 ^j. 
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has just been reached in experiment [3]. For high-Z ele- 
ments, the correction becomes large, e.g. on the level of 
10^^ for 238u9i+^ definitely will be important for the 
comparison between theory and experiment for high-Z el- 
ements, which are expected within a few years [B]. These 
results are likely to allow in future for the extraction of 
nuclear deformation parameters from experimental val- 
ues of the g factor. 

JZ and ZH acknowledge helpful conversation with 
A. I. Milstein. The work of ZH has been supported 
by the Alliance Program of the Helmholtz Association 
(HA216/EMMI). 



FIG. 1: Various contributions to the g factor; &gvL stands 
for one-loop all-order QED corrections, (5(/fs for the finite- 
size eflfect, 5(jrec for the recoil correction, (5(;nd for the nuclear 
deformation correction and &g-2i\{Zaf'\ represents two-loop 
QED corrections at order (^q)" (see [29]'). 



238U, 234u, 150^^^ i42p^(j^ loogr^ and ^^Sr, respectively. 
The accuracy of the results for the ND and QED-ND cor- 
rections is such that the comparison with experimental 
values of the g factor could serve as a method for the de- 
termination of /32 and fi^ for heavy nuclei, provided that 
all other relevant corrections as well as remaining nuclear 
parameters are known with sufficient accuracy. 

Table (|ll]) contains the values of the isotopic shift of 
iJgND for various chosen isotopes. For a comparison with 
experimental results, values of other isotope-dependent 
effects are required. The leading FS effect as well as the 
mixed FS-VP effect can be found in [12], the mixed FS- 
SE correction in [3S] , the recoil effects in , and the 

NP effect in [9]. 



TABLE II: The isotopic shifts of various corrections to the g 
factor. The following notation for an isotopic shift has been 
introduced: Apx = 5s'x(^i) — 5sx(^2), where A\ and Ai are 
the first and second mass numbers, respectively, given in the 
second column, and X G {ND,ND-QED}. 



Z Isotopes ApND 



ND-QED 



38 "'"'""Sr 

150,142j^j 



60 
92 



-9.9(2.8) ■ 10"'" 1.4(3) • 10"" 
-1.50(54) • 10"** 1.1(4) • 10"'" 
-1.6(3.9) • 10"^ 1(4) • lO"'' 



As: 



ND, Total 



-9.8(2.8) • 10"'° 
-1.49(54) • 10"* 
-1.6(3.9) ■ 10"'' 



For comparison. Fig. [T] presents a plot of 8g^\) and 
some other contributions to the g factor. It is appar- 
ent that J^ND grows quickly with Z and that it becomes 
prominent for heavy ions. 

In summary, we have developed the theory of the nu- 
clear deformation correction to the bound electron g fac- 
tor and calculated this term for a wide range of elements. 
It turns out that already for ^™Sr'^ the relative nuclear 
shape contribution exceeds 10~^°, which is the level that 
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